Statistics for Poisson models of overlapping spheres 



Daniel Hug, Giinter Last, Zbynek Pawlas and Wolfgang Weil 
January 9, 2013 

Abstract 

The paper considers the stationary Poisson Boolean model with spherical grains and proposes a fam- 
ily of nonparametric estimators for the radius distribution. These estimators are based on observed 
distances and radii, weighted in an appropriate way. They are ratio-unbiased and asymptotically 
consistent for growing observation window. It is shown that the asymptotic variance exists and is 
' given by a fairly explicit integral expression. Asymptotic normality is established under a suitable 
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integrability assumption on the weight function. The paper also provides a short discussion of related 



' estimators as well as a simulation study. 

00 

■ 1 Introduction 

(~| ' We consider a stationary random closed set Z in (d > 2) which is given as a union of random balls 

of the form 

Z ■.= [jB{U,Rn), (1.1) 

n>l 



where B{x, r) is a closed Euclidean ball with radius r > centered at x e M , $ := : n > 1} is 
^ ' a stationary Poisson point process on W^, and the sequence {Rn)n>i is independent of $ and is formed 

. by independent non-negative random variables with common distribution G. Let i? be a generic random 

'sj- \ variable with distribution G. We will always assume that it has a finite 2d-th moment, that is. 



Ei?2c« < oo. (1.2) 



Definition (11.11 ) provides an important model in stochastic geometry with numerous applications in 
physics and materials science, for instance. The set Z is called a stationary Boolean model with spherical 
grains. A simulated realization for d = 2 is shown in Figure [T] 

It is a fundamental statistical problem to retrieve information on G based on an observation of Z in a 
bounded window W. Our aim in this paper is to propose and study a family of nonparametric estimators 
of G. The nonparametric estimation of the radius distribution G has been studied before; see [3 Chapter 
5.6], lfT4l or |[T6l . In lfT6l a kernel estimator is obtained by the method of tangent points. The asymptotic 
properties of this estimator are studied in For earlier work on statistics for the Boolean model, we 
refer to [19 Chapter 3.3], [15] and the references therein. 

In the following, we assume that all random elements are defined on an underlying probability space 
($7, J", P). For a Borel set A C W^, we write ^{A) := card{n > I : Cn & A} and assume that $ has a 
positive and finite intensity 

7:=E<I>([0, 1]'^). 



MSC 2010 subject classifications. Primary 60D05, 60G57, 52A21; secondary 60G55, 52A22, 52A20, 53C65, 46B20, 
62G05. 
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Throughout the paper, let B he a compact convex set which contains the origin o and a non-degenerate 
segment. We call B structuring element or gauge body, but we point out that B need not be centrally 
symmetric or full-dimensional. The i?-distance from a point x G M'^ to a set A C M'^ is 

dsix, A) := mf{r > : (x + rS) n vl / 0} G [0, oo]. 

Clearly, if o G int B and A is nonempty and closed, then the infimum is a minimum. The most common 
structuring element is the unit ball B{o, 1), for which we also write B'^ and which is based on a scalar 
product and a norm denoted by || • ||. For given x ^ Z, almost surely dB{x,Z) < oo whenever R 
satisfies F{R > 0) > 0. We always assume that this condition is fulfilled. Then almost surely there is a 
unique n G N (that is, a ball Rn)) such that {x + ds^x, Z)B) n -B(Cn, Rn) / (see HOl Lemma 

3.1] or EOl Lemma 9.5.1]). In this case, we define rs(x, Z) as Figure [T] illustrates the definition of 
dsix, Z) and rsix, Z) fox x ^ Z and B = B'^. 




Figure 1: A simulated realization of a planar stationary Boolean model Z with spherical grains observed 
in a rectangular observation window. The symbol d denotes (x, Z) and r stands for (x, Z). 

For s, r > 0, we write B^^r '■= sB © rB'^ for the Minkowski sum of sB and rB'^. Let \A\a denote 
the (f-dimensional Lebesgue measure of a set yl C R'^, let Kk := \B''\k = tt''/^/T{1 + k/2) denote 
the volume of the /c-dimensional unit ball and write Vj{B) for the j-th intrinsic volume of B (see Il20l 
Chapter 14.3]). Then, for t G M+ := [0, oo), the empty space function Fb of Z is given by 

FB{t) ■.= ¥{dB{o,Z)<t) =P(ZntB/0) = l-exp{-7E|5t,^jU} 

= 1 - exp |-7^ | . (1.3) 

The empty space function is a useful summary statistics of random sets (see |fT9l |5l). In the case of 
a strictly convex gauge body B a detailed study of Fb for (non-stationary) germ-grain models can be 
found in ||9j|. We denote the complementary empty space function by FB{t) := 1 — FB{t). As a 
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consequence of ifTOl Theorem 3.2], we get for all measurable functions g : [0, cxd] x R+ — R"*" such that 

g(0, r) = g{oo, r) = 0, r G R+, and all x that 

/"OO /"OO 

Eg{dB{x,Z),rB{x,Z)) =j / g{t,r)hB{t,r)FB{t) dtG{dr) (1.4) 

Jo Jo 

with 

d-l 

hB{t,r) := ^(i + l)Kd-i-jVj+i{By-'-H^ , 

3=0 

for r € [0, oo); see also f20', Theorem 9.5.2]. Note that on the left-hand side of (11.41 ) the restriction to 

{0 < dB{x, Z) < oo} is expressed by the condition ^(0, r) = 5(00, r) = 0. 

For Borel sets C c M+, A C R'^ and a measurable function / : [0, 00] R+ with /(O) = /(oo) = 
0, we define a random measure r]A by 

7]a{C):= [ l{rB{x,Z)eC}f{dBix,Z))hB{dBix,Z),rB{x,Z)y^dx, (1.5) 
J A 

where !{•} denotes the indicator function. Here we put 0/0 := 0. Thus, in particular, the integration 
effectively extends over the complement Z^ := {x ^ R.'^ : dB{x, Z) > 0} of Z. Throughout the paper 
we shall assume that 

roc 

0</3:= / f{t)FB{t)dt< 00. (1.6) 



Jo 

In view of (11.31) this is a rather weak assumption on /. Moreover, we assume that the origin is an interior 
point of B if P(i? = 0) > 0. This assumption ensures that hB{t, r) > for t G (0, 00) and G-almost all 
r G R+. By Fubini's theorem and (II. 41 ). we obtain 

Er?A(C) =7/3|^|rfG(C). (1.7) 

Consider a compact convex observation window 1^ C R'^ with \ W\d > 0. We propose an estimator 
G for G based on the information contained in the data { [dB{x, Z),rB{x, Z)) : x ^ W \ Z|. Note that 
these data may also require information from outside W. The estimator is given by 



where C C R^ is a Borel set. If the denominator in (II.8I 1 is zero, then the numerator is zero as well, and 
we use the convention 0/0^= 0. From (O) we see that Er]w{C) = 7 /3 G(C7) and Erjwi'^'^) = 
1 P\W\d- This means that G is a ratio-unbiased estimator of G. 

The paper is organized as follows. In Section |2] we study second order properties of (II. 5t . We show 
that the asymptotic variance exists and is given by a fairly explicit integral expression. Consequently, 
the estimator (11.81 ) is asymptotically weakly consistent as the compact convex observation window W is 
expanding. Strong consistency follows from the spatial ergodic theorem. Section |3] contains the proof of 
asymptotic normality under an integrability assumption on the function /. In Section |4] we consider the 
estimator G in the plane and for the spherical (B = B{o, 1)) as well as for the linear case {B a segment). 
We also discuss some related estimators. A simulation study is performed to compare the behaviour of 
different discrete versions of these estimators of the radius distribution G. 
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2 Second order properties 

For a Borel set A C M°', we define the restricted Boolean model as 

Z{A):= \J B{^n,Rn). 

Clearly, Z{A) is not stationary unless A = W^. Further, for t G the complementary empty space 
function of Z{A) with respect to x € is defined by 

F^{x; t) := F(dB {x, Z{A)) > t) = P((x + tB) n Z{A) = 0) 



= E (1 - l{{x + tB) n B{U,Rn) / 0}l{Cn G A}) 

n>l 

= exp|-7E / l{{x + tB)nB{y,R) j^(l}}l{y e A}dy} 

= exp{-^E\{x + Bt,R)nA\J . (2.1) 

In particular, we have F^''{x; t) = Fsit). 

For Borel sets Ai,A2 C M'^ and ti,t2 G M+, it will be convenient to introduce the complementary 
second-order empty space Junction with respect to xi,X2 G K"' as 

F^'''^%xi,X2;h,t2):=F(dB{xi,Z{Ai)) > fi, (xs, ^(^2)) > ^2) (2.2) 



= + hB) n = 0, (X2 + ^2^) n Z(A2) = 0) 

= E J] (1 - + hB) n BiU, Rn) / 0}l{^n G Al}) 

n>l 

X (1 - 1{(X2 + t25) n Rn) + 0}l{en G ^2}) 

= exp {-7E| ((xi + Bt^^B) n Ax) U ((x2 + Sf^,/?) n A2) j^} . (2.3) 

This function is related to the second-order contact distribution function which is studied in fll. 

In order to obtain a more concise statement in the subsequent Lemma [ZT] (and again in the proof of 
Theorem l3.1l ). we introduce for given Borel sets Ai,A2 CW^ two functions, Ii (^1 , ^2 ) and I2 {Ai , ^2 ) > 



depending on the arguments {xi,X2, y,r) £ (M'^)^ x M+ and {xi,X2,yi,y2,ri,r2) G (M'^) 
respectively, which ai^e defined by 

h{Ai,A2){xi,X2,y,r) := l{y G Ai n A2}F^''^^ (^xi,X2;dB{xi, B{y,r)) ,dB{x2, B{y, 

and 

l2{Ai,A2){xi,X2,yi,y2,ri,r2) := l{yi G Ai}l{y2 G A2} 

\-l{y2 e Ai}l{dB{xi,B{y2,r2)) < (xi, n)) } 



X 



1 - l{yi G A2}1 {dB{x2,B{yi,ri)) < dB{x2, B{y2,r2) 

X F^''^^ (^xi,X2; dB {xi,B{yi,ri)) , dB {x2, B{y2, r2)) 

- (xi; (is (xi, ri)))F^2 (x2; dB (x2, B{y2, r2) 
If the arguments of these two functions are clear from the context, they are sometimes omitted. 
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Lemma 2.1. Let ^1,^2 C be Borel sets and let xi,X2 S W^. If g : [0, 00] x M+ R+ is a 
measurable function with g{0,r) = g{oo,r) = Oforr G M^, then 

Cov(^g{dBixi,Z{Ai)),rBixi,Z{Ai))),g{dBix2,ZiA2)),rBix2,ZiA2)) 

= 7 / / 9{dBixi,B{y,r)),r)g[dBix2,B{y,r)),r)li{Ai,A2)ixi,X2,y,r)dyG{dr) 
Jo JR-i 



/ / g{dB{xi,B{yi,ri)),ri)g{dB{x2,B{y2,r2)),r2) 

X l2iAi,A2){xi,X2,yi,y2,ri,r2) dyi dy2 G(dri) G(dr2). 
Proof Forn G N, 2; G W^, and i G {1, 2}, we define the event 

D^:\x) := {dB{x, \J B{^k,Rk))>dB{x,B{^n,Rn))}. 



Then 



Eg{dB{xi,Z{Ai)),rB{xi,Z{Ai)))g{dB{x2,ZiA2)),rB{x2,Z{A2))) 



= ^ X] -'-I?<,''(xi)nDi2'(x'2)^(^'^(^l' Rn)),Rn)g{dB{x2,B{^n, Rn)),Rn) 

eAinyl2 

■'■Di'>(zi)nD<^'(x2)^(^^*^^l' Rn)),Rn)g{dB{x2, B{^rn, Rm)), Rm) 

= : Ji + J2. 

Applying Mecke's formula (see [20, Corollary 3.2.3]), we obtain 



Ji=l / g{dBixi,B{y,r)),r)g{dBix2,B{y,r)),r) 
JO JAinA2 

X ¥(dB{xi,Z{Ai)) > dB{xi,B{y,r)),dB{x2,ZiA2)) > dB{x2, Biy,r))^ dyG(dr) 

and 

-^2=7^/ / / / 9{dBixi,B{yi,ri)),ri)g{dBix2,B{y2,r2)),r2) 



10 Jo JA2 JAi 

X El {dB{xi,Zy^{Ai)) > dB{xi,B{yi,ri))}l {ds (x2, Zy, (^2)) > dB{x2, B{y2,r2))} 
X dyidy2G(dri)G(dr2), 
where ^^2(^1) = ^(^1) U ^(^2,^2) if y2 G ^1 and ^^2(^1) = Z{Ai) if y2 ^ ^i- Analogously, 



5 



ZyM2) = Z[A2) U B{yi,n) if yi G A2 and ^^1(^2) = ^(^2) if yi ^ ^2- Hence, 

f'OO roc r r 

^2=7^/ / / / 9idBixi,B{yi,ri)),ri)g{dB{x2,B{y2,r2)),r2) 
Jo Jo J A2 JAi 

X (^1 - l{y2 e Ai}l{dB{xi,B{y2,r2)) < dsixi, B{yi,ri 

1 - l{yi G A2}l{dB{x2,B{yi,ri)) < dB{x2, B{y2,r2 

xF(^dB{xi,Z{Ai)) > dB{xi,B{yi,ri)),dB{x2,Z{A2)) > dB{x2, B{y2,r2) 
X dyidy2G(dri)G(dr2). 

Finally, 

Eg{dB{xi,Z{Ai)),rB{xi,Z{Ai))) = E ^ l^^^^^^^j{dB{xi, B{^n, Rn)), Rn) 

= 7/ / 5(2/1, n)), dB(2;i, ri))) dyiG(dri). 

Jo JAi 



□ 



Our aim is to analyze the second-order properties of the random measure rjA given by (I1.5I ). For this 
reason, we work with the complementary second-order empty space function (I2.2I ). For Ai = A2 = M'', 
ti, t2 G 1^^> and u = 2:2 — xi, by the stationarity of Z this function turns into 

Fj^\u;ti,t2) :=F{dB{o,Z) >ti,dB{u,Z) > ^2) 

= exp|-7E|5t,,R U (n + J 

= Fb (ii ) Fb {t2) exp {7EKB (-u; ti , t2 , i?) } , (2.4) 

where 

KB{u;ti,t2,r) := |-Bj^,rn (n + Sjj^r)]^- (2.5) 
Observe that for any u G and ti, t2 G 1^^, we have 

Fg\u;ti,t2) > FB{ti)FB{t2) (2.6) 

and 

F^b\u-MM) < exp[-^K{\Bt,,RU + \Bt,,RU)} = /S^O^^M*^- (2.7) 
These inequalities will be used subsequently. In addition, we shall need the assumption 



/•oo 

/ /(t)e-^*dt < 00, 
Jo 



(2.8) 



where c := 4 ^7Krf_iyi(i?)Ei?^ ^ < cxd and c > since Vi(i3) > (recall that B contains a non- 
degenerate line segment) and F{R > 0) > 0. 
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Proposition 2.2. Assume that (12.81 ) is satisfied. If C C M.'^ is a Borel set and Wi,W2 C are compact 
convex sets, then 



Coy {rjwAC), Vw.iC)) = / n - u)|47Ti(C, n) + tSC^, n)] du, 

where 

Ti{Cu)-= I I -^(Q' ^))) fi^Bju + X, B{o, r))) ^^^^ 

Jc Jr'' hB{dB{x,B{o,r)),r) hB{dB{u + x,B{o,r)),r) 



X F^^\u;dB{x,B{o,r)),dB{u + x,B{o,r))) dxG(dr) 



and 



r2iCu)-=f III f{dB{xi,B{o,n))) f{dBix2,B{o,r2))) 

Jc Jc Jr'' Jr-^ hB{dB{xi,B{o,ri)),ri) hB{dB{x2,B{o,r2)),r2) 

X q{u; xi,X2,ri,r2) dxi dx2 G(dri) G(dr2), 

for u G W^, and 
q{u;xi,X2,ri,r2) 

:= l{dB{x2,B{u,r2)) > dB^xi, B{o,ri))]l{dB{xi, B{-u,ri)) > dBix2, B{o,r2))] 

X F^^ {u; dB{xi,B{o, ri)),dB{x2, B{o, r2))) - Fb (dB{xi,B(p, ri)))FB{dB{x2, B{o, r2))) , 
for xi,X2 G K'^ and ri,r2 G M^". 

Proof. To abbreviate the notation, we define the function 

g{t, r) := l{r G C}f{t)hB{t, r)-\ (2.11) 

for t G [0, 00] and r G M+, with the previous conventions in the cases where t G {0, 00}. Recall also 
that hB{t, r) > for t G (0, 00) and G-almost all r G M"*". Using Fubini's theorem and stationarity, we 
get 

CoY{r]wAC),vw2{C)) 

= / Cov(g{dB{xi,Z),rB{xi,Z)),g{dB{x2,Z),rB{x2,Z)))dx2dxi 

JWl JW2 ^ ' 

= I \Win{W2-u)UCov(g{dB{o,Z),rB{o,Z)),g{dB{u,Z),rB{u,Z)))du. 

jRd ^ ' 

By Lemma [2?T] with Ai = A2 = W^, xi = o and X2 = u, we obtain that 

Cov(^g{dB{o, Z),rB{o, Z)) ,g{dBiu, Z),rB{u, Z))^ = Ji(n) + J2i(ti) - J22, 

where 

Ji{u):='y / g{dBio,B{x,r)),r)g{dBiu,B{x,r)),r) 
Jo Jr'^ 

X F^\u; dBio, B{x,r)),dB{u, B{x,r))) dxG(dr), 
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oo roo 



J2i(n):=7 / / / / 9idB{o,B{xi,ri)),ri)g{dB{o,B{x2,r2)),r2) 
Jo Jo Jr'' Jr'' 

X l{dBi-u,B{x2,r2)) > dBio,B{xi,ri))} 
X l{dB{u,B{xi,ri)) > dBio,B{x2,r2))} 

(2), 



X 



F)j {u; dsio, B{xi,ri)),dB{o, B{x2,r2))) dxi dx2 G(dri) G(dr2), 



and 



POO PCO P P 

^22:=7^/ / / / 9{dB{o,B{xi,ri)),ri)g{dBio,B{x2,r2)),r2) 
Jo Jo JR-i JR'' 

X Fb {dB{o, B{xi,ri)))FB{dB{o, B{x2,r2))) dxi dx2 G(dri) ©(dra). 

Using dB{u, B{x,r)) = dB{u — x,B{o,r)) and the reflection invariance of Lebesgue measure, we 
deduce that 

POO P 

Ji{u) = ll / g{dB{x,B{o,r)),r)g{dBiu + x,B{o,r)),r) 
Jo 

X F^\u;dB{x,B{o,r)),dB{u + x,B{o,r))) dxG(dr), 

/"OO /"OO r r 

^2i(n)=7^/ / / / 9{dBixi,B{o,ri)),ri)g{dBix2,B{o,r2)),r2) 
Jo Jo Jr'' Jr-^ 

X l{dB{x2,B{u,r2)) > dB{xi,B{o,ri))} 
X l{dBixi,B{-u,ri)) > dBix2,B{o,r2))} 

X F^^^ [u; dB{xi,B{o, ri)),dB{x2, B{o, r2))) dxi dx2 G(dri) G(dr2), 

and 

noo P P 
/ / g{dB{xi,B{o,ri)),ri)g{dB{x2,B{o,r2)),r2) 
„ jRd JlRd 

X Fb {dB{xi,B{o, ri)))FB{dB{x2, B{o, r2))) dxi dx2 G(dri) G(dr2). 

The assertion now follows by recalling (12.1 It . The integrability of ti(C, •) and T2(C, •), which is ex- 
plicitly stated in (12.151 ). will be shown in the proof of Theorem 12.41 and is implied by the assumption 
CT. □ 



Remark 2.3. Recall that || • || denotes the Euclidean norm on R*^. If B = B'^ is the unit ball, then 

dBd[u, B{x,r)) = {\\x — u\\ — r)+. 



d-l 

\d-l 



j=0 

and 



= dKS + rf-\ 



KBd{u;ti,t2,r) = \B{o,ti + r) riB{u,t2 +r)\^. 
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Hence, ti{C,u) and T2{C,u) from Proposition 12.21 may be slightly simplified. In particular, then we 
have 



T2iC,u) 



/((lki||-n)+) /((||x2||-r2)+) 
Ic JcJR'iJR'i hBd{{\\xi\\ -ri) + ,ri) /iBd((||x2|| -r2)+,r2) 

1 {(11x2 - u\\ - r2)+ > (llxill - ri)+} 1 {(||xi + u\\ - ri)+ > (||x2|| - r2)+} 



xF^^}{u;{\\xi\\-n)- 
X dxidx2G(dri)G(dr2) 



F2| 



r2) 



\xi\\-n)+)Fj,,{{\\x2\ 



r2) 




fisi 



-(si + n) 



d-l 



fis2) 



{S2 + r2 



icJcJo JS'i-^ Jo Js'i-^ h^d{si,ri) hQd{s2,r2) 

^ {{\\{S2 + r2)v2 - u\\ - r2)+ > Si} 1 {{\\{si + ri)vi + u\\ - ri)+ > S2} 



X fII}{u;si,S2) - Fb4si)Fb4s2) n''-'{dvi) dsi'H''-\dv2) ds2G{dri) G{dr2 



id-i, 



d-l, 




C JC Jo 



fisi)f{s2) 



^W^-^ {dBjo, S2 + r2) n B{u, si + r2)") 
■H'i-4dB{o,S2 + r2)) 



W^-^ {dBjo, si + ri) n S2 + ri] 

^ 7^^^-1(95(0,51 + ri)) 

X dsids2G(dri)G(dr2), 



- (2) 

-?^Sd (^i;Sl,S2) - F^d(si)FBd(s2) 



where S'^^i is the unit sphere in M'^, 'H'^^^ is the {d — 1) -dimensional Hausdorff measure, and dB{x, r) 
is the boundary of B{x, r). We used that /((||x|| — r)+) is non-zero only if ||x|| > r. Then x = {s + r)v 
for s > and v G S'^-i. 

Next we state a special case of lITOl Theorem 2.1 and Remark 3.1] in the form needed in the present 
context. Let 5 : M'^ — [0, 00] be measurable, and let K,B <ZW^ be convex bodies such that o G i? and 

B are in general relative position. Since in our application, we shall only need the case K = rB'^, for 
r G M"*", the assumption of general relative position will be satisfied for any choice of B. Then we have 



/ 1{0 < dB{z, rB'^) < oo}g{z) dz 



E 

j=0 



d- 1 
j 



t^~^-o~g(z + tb) Qj.d-j [rB''- B* ■ d(z, h)) dt, 



where B* := —B and the mixed support measures Qj-d-j{fB'^; B*;-),j G {0, . . . ,d — 1}, are finite 
Borel measures on M^'^. Using [20. (14.18)] (cf. [18. (4.2.26) and (5.3.8)]) and [12. middle of p. 327], 
we obtain for the total measures Qj.d-jirB'^;B*;R'^'^) = d{fjnjVd-j{B). In particular, this yields 

for any measurable function / : [0, 00] — [0, 00] with /(O) = /(oo) = that 



f{dB{z,rB''))d2 



hB{t,r)f{t) dt. 



(2.12) 



We now turn to the asymptotic properties of the ratio-unbiased estimator (11.81 ). Our setting is similar 
to ifTSil . where all limit theorems refer to a growing observation window in M"^. More formally, we 
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Theorem 2.4. Assume that (12.81 ) is fulfilled. Then 

^ a^{C). (2.13) 



consider a sequence (PFn)neN of compact, convex sets Wn C ffi*^ such that Wn C Wn+i for all n S N 
and the inradius of Wn tends to oo as n — oo. 



r/ie asymptotic variance is finite and given by 

^^iC)=7 [ Ti{C,u)du + ^^ [ T2{C,u)du, (2.14) 

where ti{C,u) and T2{C,u) are defined in ( 12.91 ) and (12.101 ). respectively. Moreover, ifO< G{C) < 1, 
then 0-2 (C) > 0. 

Proof. Suppose that x G Wn and u G M'^. If x + u ^ Wn, then d^dix, dWn) < Hence, we obtain 

It^nld - \{X G t^n : dB4x,dWn) < \\u\\}U < Wn H {Wn " u)\d < \Wn\d. 



Thus, ifTTl Lemma 10.15 (ii)] implies that 

\Wnn{Wn-u)U 



1 for any u G 



\Wn\d 

Therefore Lebesgue's dominated convergence theorem and Proposition 12.21 yield (12.131) provided that 

/ Ti{C , u) du < oo and / \t2{C , u)\ du < oo . (2.15) 

Using (12.71 ) we have 

rr^ ^ , , , f{dB{x,rB'')) f{dBiy,rB^)) 

Ti[C,u)du < ' ' ' 




'c 



c 



hB{dB{x,rBd),r) hB{dB{y,rB<i),r) 
Fb {dB{x, rB'^))FB {dB{y, rB'^)) dx dy G(dr 
f{dB{x,rB'')) 



hB{dB{x,rBd),', 



FB{dB{x,rB'i))dx] G(dr). 



An application of (12.121) shows that 

f{dB{x,rB'')) 



and thus we obtain 



hB[dB{x,rB'^),r) V ^ Jo ^ 

Ti{C,u)du< [ ( r f{t)J¥Mt)dt] G{dr) 
Jc \Jo ^ J 

<G{C)(^J^ /(t)e~2"*dt^ <oo, 



where we use that Fsit) < exp{— 4ct} and assumption (12. 81) . 
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In order to show that 

/ \t2{C, u)\du < oo, 

we first rewrite q{u; xi, X2,ri,r2) as the diiJerence of two non-negative terms, that is, q = qi — q2 with 

qi{u;xi,X2,ri,r2) := F^\u;dB{xi, B{o,ri)), dBix2, B{o,r2))) 

- Fb {dB{xi,B{o, n))) Fb {dB{x2,B{o, ra))) , 

which is non-negative by (I2.6I ). and 

q2{u; xi,X2,ri,r2) := F^^ {u; ^^(xi, B{o, ri)),dB{x2, B{o, r2))) 

1 - l{dB{x2,B{u,r2)) > dB{xi,B{o,ri))}l{dB{xi,B{-u,ri)) > dB{x2, B{o,r2))}) , 



for u, xi,X2 G K and ri, r2 G Using (I2.4I ). (12.71 ) and the inequality 1 — e " < a, for a > 0, we get 

Qi(u;xi,X2,ri,r2) < F^^\u; dsixi, B{o,ri)),dB{x2, B{o,r2))) 

X (l - exp{-jEKBiu; dBixi, B{o,ri)),dB{x2, B{o,r2)), R)}^ 



< ^Fb {dB{xi,B{o, n)))FB {dB{x2, B{o, r2))) 

X 'jEkb (n; dB{xi,B{o, n)), dB{x2, B{o, r2)), R) . 
Moreover, the inequality 1 — (1 — a) (1 — 6) < a + 6, for a, 6 > 0, and again (12.71) imply that 



q2{u;xi,X2,ri,r2) < \/ FB{dB{xi, B{o,ri)))FB{dB{x2, B{o,r2))) 

l{dBix2,B{u,r2)) < dBixi,B{o,ri))] 

+ l{dBixi,B{-u,n)) < dBix2,B{o,r2))] 
Combining these bounds, we arrive at 

f{dBixi,B{o,n))) f{dB{x2,B{o,r2))) 



I 



\T2{C,u)\du < 



IC JC 



hB{dB{xi,B{o, ri)), ri) Hb {dB{x2,B{o, r2)),r2) 



X ^/FB{dBixl,B{o,rl)))FB{dB{x2,B{o,r2))) 
X jEkb{u; dB{xi,B{o, ri)), dB{x2,B{o, r2)), R) 
+ l{dBix2,B{u,r2)) < dB{xi,B{o,ri))] 
+ l{dB{xi,B{-u,ri)) < dBix2,Bio,r2))} 

X dxi dx2 duG(dri) G(dr2). 

The preceding expression splits naturally into three summands which will be bounded from above sepa- 
rately. For the first bound, we observe that by Fubini's theorem 



E 



/ KBiu; si, S2, R) du = E\Bs:,,R\d\Bs2,R\d- 
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Then we apply (12.121 ) to get 

f{dB{xi,B{o,n))) f{dB{x2,B{o,r2))) 



IC JC 



hs {dnixi , B{o, n)), n) hs (^5(2:2,-6(0, r2)),r2) 



X ^JFB{dB{xl,B{o,rl)))FB{dB{x2,B{o,r2))) 
X (u; dB{xi,B{o, ri)), dB{x2, B{o, r2)), R) 



X dxi dx2 duG(dri) G(dr2) 

= 7G(C)2 /(tl)^FB(tl)/(t2)^i^i?(t2)E|St,,jjJd|i?t2,R|d dti dt2. 

Choose > such that B C cbB'^. Then |-Bt,_R|d < Kd{cBt + ii)'^ and hence the Cauchy-Schwai^z 
inequality, the convexity of s 1— )• s^, p > 1, and y/a + b < ^/a + a,b > 0, yield that 



<C3(l + tf)(l + t^), 

where ci, C2, C3 denote finite constants independent of the expectation or ii, ^2- From this and (12.81 ) it 
follows again that the first summand is finite. 

Since dB{x2, B{u,r2)) < h if and only if u S X2 + Bt-^^r-r,^ applying Fubini's theorem and (12.121 ) 
(twice) we obtain for the second summand that 




'c JC 



f{dB{x,,Bio, n))) f{dBix2,B{o, ra))) 
hB{dBixi,B{o,ri)),ri) /is (dB(x2, S(o, r2)), r2) 



'Fb {dB{xuB{o, n)))FB {dB{x2,B{o, rs))) 

X l{dB{x2,B{u,r2)) < dBixi,B{o,ri))} 
X dxi dj;2 dnG(dri) G(dr2) 

^ Lie I" I" /(ii)/^MiO/(i2)/^^|i?ti,r2ld dh dt2 G(dri) G(dr2), 

which is finite by the same reasoning as above. 

The third summand can be treated in exactly the same way. 

To prove positivity of the asymptotic variance we use the fact that the variance of any square- 
integrable function of the Poisson process ^' := {(^n, Rn) : n > 1} satisfies the inequality 

POO r 
Jo JR'i 

see, e.g., fT3^. Theorem 4.2]. In our case this means that 

Yarr]w{C)>-f / %, r)^ dy G(dr), (2.16) 

Jo JM'* 
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where 

h{y,r) :=E [ [g{dBix, Z U B{y,r)),rB{x, Z U B{y,r))) - g{dB{x, Z),rB{x, Z))] dx 
Jw 

= E / l{dBix,B{y,r)) < (is(x,Z)} [g{dBix, B{y,r)),r) - g{dB{x, Z),rB{x, Z))] dx 
Jw 

= E / l{dBio,B{y-x,r)) <dBio,Z)} 
Jw 

X [9{dB{o,B{y - x,r)),r) - g{dB{o, Z),rB{o, Z))] dx. 
Here the last identity follows from the stationarity of Z and g is as defined in (12.1 lb . By (11.41 ). 



oo roo 



h{y,r)=-f / / l{dB{o,B{y-x,r)) <t} 
Jw Jo Jo 

X [g{dBio, B{y - x, r)), r) - g{t, s)] hB{t, s)FB{t) dt G{ds) dx. 

Assume now that < G(C) < 1 and let C := M+ \ C. Recalling the definition (12.111) of g, we obtain 
from (12.161 ) that 

/"OO r 

Varr/iy(C7) > 7 / / h*{y,rfl{r £ C' ,y £ W} dyG{dr), 
Jo JR'' 



where 

h*{y,r) :=7 



f' f'OO f'QO 

/ / / l{dBio,B{y-x,r)) <t}g{t,s)hB{t,s)FB{t)dtG{ds)dx 
Jw Jo Jo 

7G(C7) / / l{dB{o,B{y-x,r))<t}f{t)FB{t)dtdx. 
Jw Jo 



Applying Jensen's inequality with the normalization of l{r G C", y G M^} dy G(dr), we get 

, ^. ,,,,,, , ^ , , ... h* (ii r \ c\ii (Hri c\r\ 

'-GiC'WU KJc'Jw 
Letting a := 'y^G{Cf /G{C') > we obtain that 



\wu 



a If r I r iS^dB{o,B{y-x,r)) <t}f{t)FB{t)dtdxdyG{d ^ 
Jc Jw Jw Jo 

moo 
\WniW- y)\dl{dB{o, B{y,r)) < t]f{t)FB{t) dtdyG{d 



a 



Hence it is sufficient to show that 



/ l{dB{o,B{y,r)) <t]dy\ f{t)FB{t)dtG{dr)>{). 

Jc Jo xJw J 

This is true, since the inner integral is positive for all r, t > and since both \^ f{t)FB{t) dt and G{C') 
are positive. □ 
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Remark 2.5. The assumption (12.8b is slightly stronger than (ll.61 l. 

Remark 2.6. Let G„(C) be given by (fTSl) with = W„. Theorem IZ4] implies that G„(C) is asymp- 
totically weakly consistent. Indeed, (11.71 ) and 

n-s>oo 

ensure that r?vF„ (C*) /|^n|d converges to 7 /3 G(C) in probability as n — 00. Especially, 



n\d 



7/3 in probability . (2. 17) 



Hence, by the continuous mapping theorem, rjw^^C) /'qwi'^'^) converges to G{C) in probability as n — )• 
00. This is in accordance with the following proposition which even shows that Gn{C) is asymptotically 
strongly consistent. 

Proposition 2.7. For any Borel set C C M+, we have G„(C) — > G(C) ¥-a.s. 

Proof. The mapping W 1— )• r]w{C) defined by (11.51 ) is a random measure on M'^ depending on the 
Boolean model Z in a translation-invariant way. As the Boolean model is ergodic (see ll20l Theorem 
9.3.5]) we can apply the spatial ergodic theorem (see iflTl Corollary 10.19]) to conclude that 

lim iWnlfvwAC) = Er/[o i]d(C) = 7/3G(C) P-a.s. 
Applying this to the numerator as well as to the denominator in (11.81 ). we obtain the desired result. □ 



3 Asymptotic normality 

In this section we study the asymptotic normality of the ratio-unbiased estimator (11.81 ) for the radius 
distribution G of our stationary Boolean model Z with spherical grains. The proof will be based on 
approximation by m-dependent random fields. This idea comes from [7], where the same technique 
was used to prove the central limit theorem for random measures which are associated with the Boolean 
model in an additive way. In contrast to [7], the contribution of an individual grain to the random measure 
A I— >• r]A{C) is not determined by the grain alone, but does depend on a random number of other grains 
in a non-trivial manner. Therefore the results of [7 1 do not apply in our setting. 
We consider, forn G N and a Borel set C C M+, the estimator 

where Wn ■= [— n, n)"^ and rjWn given in (11.51 ). First we concentrate on the asymptotic normality of 
the numerator r/vi/„ (C)- In addition to (11.61) . we shall need the integrability condition 

/•oo 

/ (1 + dt < 00, (3.1) 

Jo 

which is more restrictive than (12.81 ). 

Theorem 3.1. Assume that (11.61 ) and (13.11 ) are fulfilled. Then, for any Borel set C C M"*", 

V Wn\d J 

where cr^(C) is given by (12.141 ). 
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Proof. We fix a Borel set C C K"*" and skip tlie dependence on C in the notation. Let Ez := [0, l)*^ + z 
for z € In '■= {—n, . . . ,n — 1}'^. Then 



XI = X] / 9{dB{x,Z),rBix,Z)) dx, 



Z£lri zel, 



where g is given by (I2.11I ). For some fixed integer m, we put := Ez © [—m, mY- We decompose 
TjE^ into two random variables 

r?f):= / g{dB{x,Z{Fz)),rB{x,Z{Fz)))dx 



, ~(m) (m) J . (m) v-^ (m) , ~(m) v-^ ~(m) . (m) , ~(m) 

and 7]'z := 7]e, - r]l ' . Let r/^j := ^^^^^ r/^ and r]\^^ := }2z(^i„ so that 7/iy„ = + r/^^;. 
It is easily seen that {ryi™'' : n G [/} and {r/i'"'' : v € y} are independent whenever U,V dU^ are such 
that ||u — ^lloo > 2m for each u ^ U and v G y. Thus, the random variables r]^z^\ for z G Z"^, constitute 
a stationary (2m)-dependent random field (cf. |0 Section 4.3.1]). The variance of -q^^ is 

where Nn{z) is the cardinality of {(^1,22) ^ In x In '■ Z2 — zi = z}, which may be bounded by 
\Wn\d = {"^nf and lim„^oo Nn{z)/\Wn\d = 1 for any z e U^. We define 

I l^nld 

Since Cov(r/o™\ r/i™"*) = for ||2;|| > 2m, the limit of (ai™^)^ as n — )• 00 exists and satisfies 

:= lim = J] Cov(7?(™),7?(™)). (3.2) 

n— s-oo ^ — ' 

2e{-2m,...,2m}'* 

Next we show that E(r7o"^^)^ < 00. We put A := [—m, m + 1)*^, hence 

Vir^= I g{dB{x,Z{A)),rB{x,Z{A)))dx. 

Proceeding as in the proof of Proposition I2.2l and bounding F^'^{-) as well as (1 — 1{-}1{-}) by 1, we 
get 

E(r/("'))' <7 / / / [ g{dB{y,B{xi,r)),r)g{dB{y,B{x2,r)),r)dyG{dr)dxidx2 
JEo J Eo Jo J a 



p p poo poo p p 

+ 7' / / / / / / g{dBiyi,Bixun)),n) 

J Eo J Eo JO Jo J A J A 

X g{dB{y2,B{x2,r2)),r2) dyi dy2 G(dri) G(dr2) dxidx2. 

The right-hand side increases if A is replaced by M*^. Arguing then as in the proof of Theorem 12.41 we 
obtain 

(foo \ 2 / /.oo \ 2 

f{t)dt] G{C)+ h\EoUJ^ /(t)dtG(C)j <oo. 
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Therefore, the central limit theorem for stationary m-dependent random fields (see, e.g., UTl) yields that 



In view of |2, Theorem 3.2], it remains to verify that 

lim = a{C) (3.3) 



and 



Define 



lim lim sup P ( , ^ I V (Ty^™) - Er/i™)) > e | = for any e > 0. (3.4) 



2 Var r]w„ 



\Wn\d 

Then 

W{C) - ^7(-)| < \a{C) - an\ + Wn - + ki'") - 

The first term goes to zero as n — )• oo by Theorem 12.41 the last term goes to zero as n — )• oo as well, for 
any m G N, by (13.21 ). By Minkowski's inequality, the middle term can be bounded as 



Therefore, (I3.3l i follows if we can show that 



By Chebyshev's inequality, (13.51 ) also implies (13.41) . The variance in (13.51) satisfies 

-^Var j;r/f)=-^ J] iV„(z) Cov(r?(-), r?f)) < ^ | Cov(r?(™), ^M) |. 

Therefore, the proof will be finished when we show that 

|Cov(,?(™),,7i-))| 0. 

Consider a fixed z G Z'^. Then the covariance can be written as 

Cov(r?(™),r?f)) = Coy{r^E.,r^E.) - Gov (,?(™) , ,7^ J - Gov , r?(-) ) + Cov(ry(-) , r?^)) 

= / / [cRd,iRd(a;i, 3^2) - c^„,R<i (a:i,X2) - Cigd^^^(xi,X2) + CF„,F,(2;i,X2)] dx2dxi, 

where, for Borel sets Ai,A2 C M'^ and xi, X2 € M*^, 

CAi,A2(2;i, X2) := Cov(5(dB(xi, Z(^i)), rB(xi, .^(^1))) ,g{dB{x2, Z{A2)),rB{x2, -^(^2))) 
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is expressed in Lemma llTI as 



oo 

Jw^ 

oo roc 



g{dB{xi,B{y,r)),r)g{dBix2,B{y,r)),r)li{Ai,A2)dyG{dr) 
g{dB{xi,B{yi,ri)),ri)g{dB{x2, B(y2, ^2)), 



CAi,A2ixi,X2) = 7 



7 

JO JO 

X 12(^1,^2) dyi dy2G(dri)G(dr2). 

Here we skip the arguments xi,X2, y,r, respectively xi,X2,yi,y2,ri,r2, of tiie functions Ii{Ai,A2) 
and 12(^1,^2), wiiich were defined before Lemma ITTI We shall treat both parts of cai.AjC^I' ^2) 
separately. Our aim is to prove that 



and 



'Eo JE^ Jo 

ri(M"',M'^) 

00 POO 



a-EoJE,JO Jo 
d -md 



h 



h{Fo 



g{dB{xi,B{y, r)),r)g{dB{x2, B{y, r)),r) 
Ii(Fo,M^) -Ii(M^F,) +/i(Fo,F,)| dyG(dr)dx2dxi 

g{dB{xi,B{yi,ri)),ri)g{dB{x2,B{y2,r2)),r2) 
^) - /2(M^ F,) + hiFo, F,) dyi dy2 G(dri) G(dr2) dx2 dxi 



tend to zero as m — )• 00. Observe that Si, S2 depend on m via the dependence of Fq, Fz on m. 
First, we consider Si . We rewrite 



h( 



hd Tnid\ 



hiFo, M'^) - Fz) + /i(Fo, Fz 



l{y G F^ n F,^}F^"'"-"(xi, X2; ti, ts) + Hv G Fo n F,^} ( F^"'"" - F^°'"" ) (xi, X2; ti, t2) 



+ l{yGF„=nFj^F^ 
+ l{y G Fo n FJ (f^ 



F 



{xi,X2;ti,t2) 



(xi,X2;tl,t2) 



with ti = dB{xi, B{y,r)) and ^2 = dB{x2, B{y,r)y For notational simplicity, write 

z^i(xi,ti) := jEUxi + Bt,,R) n F^\,, V2{x2M) := 7^1 (3^2 + ^t^,/?) n F 



^ Id 



for xi,X2 G M*^ and ti,t2 £ l^^- We suppress the dependence on z in i'2{x2,t2)- From (12.31 ) and the 
inequality 1 — e~" < a, for a > 0, we obtain for xi,X2 G M'^ and ti, ^2 G that 



F 



{xi,X2;ti,t2) 



Fi 



(xi,X2;tl,t2) 



= F^°'^'\xi,X2;ti,t2) (1 - exp {-7E| (xi + 5*,,^) n (x2 + Bt,,RY n F„^|J 
< ui{xi,ti). 



Analogously, 



Fg ' - Fg ' M (xi,X2;tl,t2) <J^2{x2,t2)- 



(3.6) 
(3.7) 
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Furthermore, 



B 



1 - exp {-7E|(xi + Bt,,R) n n ((x2 + Bt^^Rf U 
1 - exp {-7E|(X2 + Bt,,B) n n ((xi + 5^,,^)^ U F^)\j') 
+ 3^2; ti, i2) (1 - exp {-7E| (xi + Bt,,R) D (x2 + ^^.r) n D F^\J^ (3.8) 



gives 



< iyi{xi,ti)iy2ix2,t2) + \/z^i(xi,ti)i^2(a;i,i2) 
because by the Cauchy-Schwarz inequaUty 



(3.9) 



E|Xi n X2\d < Vnx^dVnx2\d 

for any random sets Xi and X2. Combining (13.61) . (13.71 ) and (13.91) . we obtain 



Tnia!\ 



hiFo, 



/i(M^F,) + /i(Fo,F, 



< l{y G n F,^} + l{y G Fo n F,^}i/i(xi, ti) + l{y G n Fjzy2(x2, t2) 



+ l{yGFonF2} I^l{xi,ti)l^2ix2,t2) + y^i^iixi,ti)i^2{x2,t2) 



(3.10) 



where ti = dB{xi, B{y,r)) and ^2 = dB{x2, B{y,r)y If xi G Eo and X2 G F^, then i/i{xi,ti) is 
bounded by i^(ti) and ^'2(x2, ^2) is bounded by z/(t2), where 

:=7E|(Foe5t,fi) nF„^|^, tGM+. 

Let Cfi > be such that S C c^^"'. Then 

zv(t) <E|5^_^^^|^<KdE(cBt + i? + ^/d)'^<ci(l + t^), (3.11) 

where ci is a finite constant that does not depend on t. If xi G Fq and y G F^, then ||xi — y\\ > m 
and dB{xi, B{y,r)) > [m — r)'^ /cb- Similarly, if X2 G F^ and y G F^, then dB{x2, B{y,r)) > 
(m — r)~^ /cb- Let 

V'(t,r) := 1 <^ t > ^ ^ 



, t,r G 



Then, by (13.101 ) and the substitutions xi — y — )• xi and X2 — y — ?■ X2, we get 

Si<y\ [ [ [ [ g{dB{xi,rB'^),r)g{dB(x2,rB'^),r)l{xi + yeEa,X2 + yeE,} 

-'K'' •^K'* JO JR'i 

^idB{xi,rB'^),r)^lj{dB{x2,rB'^),r) + V((iB(x2, rS'^), r)z.(dB(xi, rS'^)) 
+ ij{dB{xurB'^),r)i^{dB{x2, rB'^)) + u{dB{xi,rB'^))v[dB{x2,rB'^)) 



X 



+ \/i'{dB{xi,rB'^))i'[dB{x2,rB'^)) dy G(dr) dx2 dxi, 
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and thus 



Si < 



JR'i 



g{dB{xi,rB'^),r)g{dB{x2,rB'^),r) 



X V 



{dsixurB'^), dB{x2, rB'^),r) dx2 dxi G(dr 



where 



v{h,t2,r) := {iPiti,r) + i^{ti)) (V'(t2, + zv(t2)) + VHhHt2), h,t2,r G M+. 
Now Fubini's theorem, two appUcations of (12.121 ) and definition (12.111 ) of g yield 

Si< I I I fih)f{t2)v{h,t2,r)dt2dhG{dr). 



c Jo Jo 



Our moment assumption (11.21 ) ensures that E|£^o © ^t,R\d ^ °° ^'^^ therefore u{t) — as m — oo 
by Lebesgue's dominated convergence theorem. By (13.1 11 1 and (13. II ). another application of Lebesgue's 
dominated convergence theorem shows that Si — as m — oo. 

Next, we proceed with S2. From (12.41 ) and the inequality 1 — e""^ < a, for a > 0, we obtain for 

xi,X2 e W^' and ti,t2 G M+ that 



F^ '^ (xi ,X2;h,t2)- (xi ; ti (^2 ; t2) 

= '"^ {xi,X2;ti,t2)(l - exp{-'yEK.Bix2 - xi;ti,t2,R)} 
< 'y¥,K.Bix2 - xi;ti,t2,R), 
where is defined in (12.51 ). and 



(3.12) 



,„_pF.,.^^ (xi,X2;ti,t2) - [F^^-F'^°) {xi;h)F^^{x2;t2) 
= -F^°'^\xi,X2;ti,t2)(l - exp {-7E|(xi + Bt,,R) n {x2 + St^,/?,) n F^ 
+ (xi,X2;ti,t2)(l - exp {-7E|(xi + Bt.^n) n (x2 + -8*2,/?)!^ 



Therefore, 



{xi,X2;h,t2) - [F^" -F^°){xr,ti)F^^{x2;t2) 
< 7Ekb(x2 - xi;ti,t2,R). 



B 



(3.13) 



Analogously, 



B 



B 



ixuX2;h,t2)-F^ (xi;ti) F| -F^M(x2;t2) 



< 7Ekb(x2 - xi;ti,t2,R). 



(3.14) 
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Finally, using (13.81 ) we get 

-{Pr-Fi^) ixr,h){Fr-F^^)ix2;t2) 
= F^°'^' (xi, X2; ti, t2) (l - exp {-7EI (xi + Bt,,R) n F^ n ((x2 + Bt^^Rf U F,^) |J ^ 
^1 - exp {-7E|(x2 + Bt,,R) n F^ D ((xi + Bt^Rf U F^) \ 
+ F*''^'(xi,X2;ti,t2)(l - exp {-7E|(xi + Bt,,R) H (x2 + Bt,,R) H F^ n F^^ 
-F|'°(xi;ti)F^^(x2;t2)(l-exp{-z.i(xi,ti)})(l-exp{-zy2(x2,t2)}) 
= (xi, X2; ti, t2) (1 - exp {-7E| (xi + Bt,,R) D F^ D ((x2 + Bt^^Rf U F^) 

X (1 - exp {-7E|(x2 + Bt.^R) D F^ D ((xi + Bt.^Rf U F^] 

X (1 - exp {-7E|(xi + Bt.^R) D (x2 + Bt^,/?,) n Fo n F,\J^ 

+ F^''^'(xi, X2; ti, t2) (1 - exp{-7E| (xi + Bt,,R) D (x2 + Bt.^R) D F^ D F^\^} 

- F|'°(xi;ti)F^^(x2;i2) [exp {-7E|(x2 + Bt,,R} n n ((xi + Bt.^RY U F^)\J 
X (1 - exp{-7E|(xi + n (X2 + Bt^^R) n Fo n F^\^}) (l - exp{-z^i(xi, ti)}) 

+ exp {-7E|(xi + Bt^^R) n F„^ n ((x2 + St^,/?)" U F^^) [J 
X (1 - exp {-7E|(xi + Bt,,R) D (x2 + Bt^^R) D F^ n F^jJ) 

X (1 - exp {-7E|(X2 + Bt,,R) n F,^ n ((xi + Bt.^RY U F„^) |J) 
which leads to 



< Z^i(xi,ti)z^2(2;2,t2)7lE«:B(x2 - Xi]ti,t2,R) 

+ 7E|(xi + n (x2 + 5*2, r) n F^ n F^^l^ 

+ l/i(xi,ti)7EKB(x2 - Xi\ti,t2,R) + V2{x2,t2)^^KB{x2 - Xi] ti, t2, R) ■ (3.15) 

In the following, we use (13.121 ) - (13.151 ) with ti = ds(xi, ri)), t2 = dB{x2, B{y2,r2)). More- 
over, we define 

Xi := l{dB{xi,B{y2,r2)) < h} = 1 {xi G ^2 + ^ti.rj} , 
X2 ■■= 1 {dB{x2,B{yi,ri)) < 12} = 1 {x2 G yi + ^ta.n} , 
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and 

KB{xi,X2;ti,t2,R) := |(xi + Bt^^R) n (x2 + Bt^^n) n F^\^. 
Next, we fix xi £ Eq and X2 € E^, for the moment, and distinguish several cases. 

1. If 2/1 G and y2 G , then ti > (m — ri)+/cB, t2 > ("i — ^2)^/05 and using (13.121 ) we get 

/2(M^M'^)-/2(i^o,K'^)-/2(K^F,) + /2(i^o,i^.) = /2(M^M'^) 



= (1-Xi)(1-X2)Fb (2;i,X2;tl,t2)-i^B(tl)FB(t2) 
< '-f¥.KB{x2 - Xi;ti,t2,R) +X1+ X2- 

2. If yi G Fo and 2/2 G n Fq, then t2 > ("t- — '^2)^ /cb and using (13.131 ) we get 



^^>c^^ 



/2 (Fo , M"^) - /2 (M'^ , ) + /2 (Fo , ) 



^^)C^^ 



/2(Fo,M'^) 



(1-Xi)(l-X2) Fj 



(xi,X2;tl,t2) 



(xi;ti)FB(t2 



< -^'Ekb{x2 - Xi; Si, S2, i?) + Xl + X2- 

3. If yi G Fo and 2/2 G F^ n Fq, then t2> {m — r2)^ /cb and using (13.131 ) we get 



/2(Fo,M'^) - /2(M^F,) + /2(Fo,F,) 



-md\ 



h{Fo, 



(I-X2) (l-Xi)i^j 



B 



{xi,X2;ti,t2) 



B -^B 



{xi;ti)FB{t2 



< 7Ekb(x2 - xi; si, S2, -R) + xi + X2- 
4. If j/i G Fq n Fx and 1/2 G F^, then ti > (m — ri)+/cB and using (13.141 ) we get 



I2 



V^2(K^F, 



, M'^) - /2(Fo, M'^) - /2(M^ F,) + /2(Fo, F, 

= (1 - Xi)(l - X2) (f^'^' - Fr'""^) (^i,^2;ti,t2) - FB(ti) (F- - F-^ ) (X2;t2) 
< 7Ekb(x2 - xi;ti,t2,R) + X1+ X2- 
5. If yi G Fo n F^ and 7/2 G F^, then > (m — ri)+/cB and using (13.141 ) we get 



/2(Fo, M'^) - /2(M^ F,) + /2(Fo, F, 



V^2(M^F, 



(1 - Xi) (1 - X2)F^ - F^ (xi, X2; ti, t2) - FBih) Ff - F^^ (X2; t2) 



< 7Ekb(x2 - xi;ti,t2,R) + X1+ X2- 
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6. If yi E Fo n and y2 G n Fq, then by ([321) and (I3l5l) . 

/2 (M'^ , M"^) - /2 (Fo , M'^) - /2 (K^ F, ) + /2 (Fo , F, 



(1-Xi)(l-X2) 



F^'^^+F^-^^)(xi,x2;ti,t2) 



< {l^l{xi,ti)u2{x2,t2) +Z^l(xi,ti) + 1^2{x2,t2))7^KBix2 - Xi;ti,t2,R) 



+ 7EKB(xi,X2;tl,t2,-R) + (Xl + X2) i'^lixi,ti)iy2ix2,t2) + ^/'^l{xi,ti)l^2ix2,t2) 



7. If yi G Fo n and ^2 G F^ n Fo, then by ^M, ^3 and (IXBl) . 



/2(M^M'^)-/2(Fo,M'^)-/2(M^F,) + /2(Fo,F,) = (1 - - ^2) 



X F: 



ixi,X2;tut2) - (1 - Xl) (i^f"''^ - Ff°'^' 



B 



{xi,X2;ti,t2) 



< {iyi{xi,ti)iy2{x2,t2) + + Z/2(2;2,t2))7lEKB(x2 - Xi;ti,t2,-R) 



+ 'yEKB{xi,X2;ti,t2,R) + Xi {j^i{xi,ti)u2{x2,t2) + \/i^i(a;i,ti)i/2 (3:2,^2) 
+ X2i^i(a;i,ti). 

8. If yi G Fo n F^ and y2 G F^ n F^, then by (O, dlSll and dTTSl) . 



X F: 



/2(Fo,M^)-l2(M^F,) + l2(Fo,F,) = (l-xi)(l-X2) 



f;^ '^0 (xi,X2;ti,t2)-(l-X2)(FB° 



F^°'^^){xuX2;h,t2) 



Pr-F^^) {x,;h) [Pr-Fi^) {X2-M) 

< {h'l{xi,ti)u2{x2,t2) + Ul{xi,ti) + h'2{x2,t2))l^K'B{x2 " Xi;ti,t2,-R) 
+ 7E/iB(xi,X2;tl,t2,-R) + Xl'^2(2;2,i2) 
+ X2 {l^l{xi,ti)v2{x2,t2) + \/ yi{xi,ti)v2{x2,t2)^ ■ 



9. If yi G Fo n F^ and y2 G F^ n Fq, then xi = 1 implies ti > (m — r2)^ /cb and X2 = 1 imphes 
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t2 > {m — ri)+/cB, and by (13.151 ) we have 



pd jnid\ 



B 



(1 - X2)i^B°'"" - (1 - xOi^s"'"^ + Fb°'"' ) (^1' h,t2) 



< {iyi{xi,ti)iy2ix2,t2) + + iy2{x2,t2))l^KB{x2 - Xi;ti,t2,R) 

+ jERBixi,X2]ti,t2,R) + Xl + X2- 

Altogether this gives 



- /2(Fo,M'^) -/2(M^F,) + /2(Fo,F,) < /2*(xi,X2,yi,y2,r-i,r2), 



where 



l2ixi,X2,yi,y2,ri,r2) := n)V'(*2, ^^2) + 2V'(t2,?'2) +2ip{ti,ri)) 

X (7EKB(a;2 - xi;ti,t2,R) + Xi + X2) 

+ A{u{ti)u{t2) + l^ih) + u{t2))jEKB{x2 - Xi;ti,t2, R) + 4-fEkB{xi,X2;ti,t2, R) 



+ 2(xi + X2) v{ti)v{t2) + yjv{ti)u{t2) 



+ X2i^{h) + Xi'^ih) + Xi'^{h,r2) + X2^(i2, n) 

does not depend on z. We recall the dependence of XI1X2 on xi,X2, 2/1,2/2 and then carry out the 
substitutions xi — yi ^ xi and X2 — 2/2 — ^ X2 to get 

^2<E / / / / / / 9{dBixunB''),n)g{dBix2,r2B''),r2)l{xi+yiGEo} 

X l{x2 + 2/2 e -E'z}/2(2;i +2/i,a;2 + 2/2,2/i,2/2,n,f'2)d2/id2/2G(dri)G(dr2)dx2dxi 

= / / / / / / 5(dB(xi,riS'^),ri)5((iB(x2,r2B'^),r2)l{xi+2/i G^o} 
JR<* JRd JO JO JR<* JR'' 

X /^(xi + 2/i,a;2 + 2/2,2/i,2/2,n,r-2) d2/i d2/2 G(dri) G(dr2) dx2 dxi. 
To the inner integrals we apply the relations 

/ / l{xi + 2/1 G £^o}xi d2/2 d2/i = / / +2/1 - 2/2 G 5ti,r2}d2/2d2/i = l-Bti.raU, 

I I l{xi + 2/1 G £^o}X2 d2/2 d2/i = / / l{x2 + 2/2 - 2/1 G 5*2,^} d2/2 d2/i = l^t^^^i U, 

/ / l{xi + 2/1 G Eo}'Ekb{x2 + y2-xi- yi;ti,t2,R)dy2dyi = E\Bti^R\d\Bt2,R\d, 

JRd JR<^ 

/ / l{xi +2/1 G £^o}Ekb(xi + 2/i,X2 + 2/2;ii,i2,-R)d2/2d2/i 
= / l{xi + 2/1 G ^o}E|(xi + 2/1 + Bt,,R) n F^\d\Bt,,BU < E|(^o Bt^^^) n F^^UIS^^.rU. 
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In the last two equations we used |[20l Theorem 5.2.1]. Consequently, (12.121) and (12.111 ) yield 

moo /"CXD p 
/ f{h)fit2) (V'(ti,ri)V'(t2,r2) + 2V'(t2,r2) + 2^(ti,ri)) 
_ _ _ Jo L 

X {7^\Btj^,R\d\Bt2,R\d + \Bti,r2\d. + \Bt2,ri\d) 

+ 4(i^(ti)zy(t2) + iy{h) + i/(t2))7lE|5ti,ijU|5i2,RU + 47E|(So e Bt,,R) n F<$=|d|St2,ij|d 



+ \Bti^r2\dtp{h,r2) + |-Bt2,nldV'(*2,n) dt2 dti G(dr2) G(dri). 

Using E|i?f_/j|rf < EKii'offi-Bt r)!^ < ci(l + t'^), cf. (13.111 ). the Cauchy-Schwarz inequality and assump- 
tion (I3.1I ). it follows from the Lebesgue dominated convergence theorem that 5*2 — )• as m — )• 00. 

□ 

Now we are dealing with the asymptotic normality of G„(C). 
Theorem 3.2. Assume that (13.11) is satisfied. Let Wn = [— n, n)"^. If C C M+ is a Borel set, then 

^WM (g„(C) - G(C)) A iV(0,4(C)), 

where 

^UC) := ^ [(1 - G(C)y\C) + G(C)a2(M+ \ C) - G(C)(1 - G(C))a2(M+)] (3.16) 
and o-2(.) ^/ve« hy (12141) . 7/0 < G((S') < 1, f/ien cr^ (C) > 0. 



Proof. Using (12.171 ) and Slutsky's theorem, the weak limit of y^l W„|rf (^G„((7) — G((7) j coincides with 
the weak limit of 

Observing that 

7/3y„ = ^i= / {l{rB{x,Z)GC}-GiC))f{dB{x,Z))hB{dB{x,Z),rBix,Z))-Ux, 

we can proceed along the same lines as in the proof of Theorem 13. H and obtain 

7/3yn A iV (0,7^/3^2(0), 

n— 5>oo 

provided we can identify the asymptotic variance i7q(C) of y„. Theorem I2.4l implies that 

7^/3^ lim VaiYn 

n— >oo 

= a2(C) + G(C7)V(M+)-2G(C) lim ^_ Cov(wJC), w„(M+)). (3.17) 

n^oo \ Wn\d 
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Since ?7vy„(-) is additive, we obtain from Theorem l2.4l that 

2 lim-^Cov(w„(C7),w„(K+))=2a2(C) + 2 lim ^ Cov(r?w^„(C),7?w^„(R+ \ C)) 

= 2a^ (C) + lim (Var 7]w^ - Var (C) - Var r?i^„ (M+ \ C)) 

= 2a\C) + a\R+) - a^{C) - a^{R+ \ C) = a\C) + a\R+) - a\R+ \ C). 

Inserting this result into (13.171 ) we obtain (I3.16l l upon some simplification. 

To prove the last assertion, we define g{t, s) := G C} — G{C))f{t)hB{t, s)~^ and assume that 
< G(C) < 1. For a convex body 1^ C M'^ we need to consider the variance of 

Hw ■= / g{dB{x,Z),rB{x,Z)) dx. 
Jw 

As in the proof of the positivity assertion in Theorem I2.4l we obtain that 

poo r 

YaTHw>7 / h{y,rfdyGidr), (3.18) 
Jo Jr'^ 

where 

h{y,r) :=E / l{dB{o, B{y - x,r)) < dsio, Z)}g{dB{o, B{y - x,r)),r) dx 
Jw 

-E / l{dB{o,B{y - x,r)) < dB{o,Z)}g{dB{o,Z),rB{o,Z)) dx. 
Jw 

By (11.41 ) and the definition of g the second expectation on the above right-hand side vanishes for all 
y E W and r > 0. Therefore 

h{y,r)= I FB{dBio,B{y - x,r)))g{dB{o,B{y - x,r)),r)dx. 
Jw 

Again as in the proof of Theorem l2.4l we let C := \ C and obtain from Jensen's inequality and (I3.18l l 
that 



-y/Var Hw ^ c 



/ / / FB{dB{o,B{y-x,r)))f{dB{o,B{y-x,r))) 
Jc Jw Jw 



^/\W\d ~ \Wl 

X {hB{dB{o,B{y - x,r)),r))~^ dxdyG(dr) 

\Wr\{W- y)UFB {dB{o, B{y, r)))f{dB{o, B{y, r))) 



\W\ 



d Jc 



X {hB{dBio,Biy,r)),r)) 'dyG(dr), 
where c > is a constant not depending on W. Hence it is sufficient to show that 

Fb {dB{o, B{y, r)))f{dB{o, B{y, r))) {hs {dB{o, B{y, r)),r) ) dy G(dr) > 0. 
By (12.121) the above integral equals G(C') J^a (*)/(*) dt, which is positive by (11.61 ). □ 
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c 



f{dB{o,B{x,r))) f{dB{u,B{x,r))) 
Ad hB{dB{o,B{x,r)),r) hB{dB{u, B{x,r)),r) 

X F^^\u;dB{o,B{x,r)),dB{u,B{x,r))){l{r e C} - G{C))^ dx G{dr) 



Remark 3.3. After some manipulation we get 

7'/3V2(C)=7 [ hiC,u)du + ^^ I f2{C,u)du, 

where 

h{C,u) :-- 

and 

f{dB{xi,B{o,n))) f{dB{x2,B{o,r2))) 
JO JRd JR^ hB{dB{x2,B{o,ri)),ri) hB{dBix2,B{o,r2)),r2) 

X 1 {dB{x2,B{u,r2)) < dB{xi,B{o,ri))}l {ds (xi, S(-u, n)) < (x2, S(o, r2)) } 

X F^^\u;dB{xi,Bio,n)),dB{x2,Bio,r2))){l{n G C} - G(C)) (l{r2 G C} - G(C)) 
X dxidx2G(dri)G(dr2). 
This relation can be also obtained directly by an analogue of the proof of Theorem 12.41 

4 The planar case 

We mentioned at the beginning that the estimator G which we discussed so far is based on the data 
{ {dB{x, Z),rB{x, Z)) ■.x£W\Z] and therefore may require information from outside the window 
W. To overcome this problem, a common procedure in spatial statistics is the so-called minus sampling, 
which can be used, e.g., if the radius distribution G is concentrated on an interval [0, tq], < tq < oo. 
We can avoid such a condition by assuming that the function / is concentrated on an interval [0, e] with 
e > 0. If we then assume that Z is observable in a window VF^^) which contains W ® eB, then, for each 
X G W, we have either f{dB{x, Z)^ = or dB{x, Z) < e, in which case the (almost surely unique) 
contact point (x + c?b(x, Z)B) n Z lies in VF^^). 

For practical applications, the planar case d = 2 is particularly important. Also, then, the spherical 
case B = B^ and the linear case B = [0, u] (with a given direction u) play a major role. For simplicity, 
in the following considerations we concentrate on the window W = [0, 1]^ and we assume, as explained 
above, that / is concentrated on [0, e],e > 0, and that Z is observed in W^^h Let Ci, . . . , (7^ be the 
(connected and relatively open) visible arcs in dZ n VF*-^^. We need not know whether some of these 
arcs belong to the same particle. By Ci C Ci, . . . ,Ck C Ck denote the corresponding "effective" 
arcs; these consist of the points (x + dB{x, Z)B) n Z £ Ci, x £ W \ Z,for: which dB{x, Z) < e. Let 
Tj be the radius and li the length of Ci, and let Ai be the set of points x € W\Z with dB{x, Z) < e and 
which project onto Ci, i = I, . . . , k, in the sense that {x + dB{x, Z)B) n Z consists of a unique point 
and this point lies in Ci, for x & Ai. Then our estimator G is of the form 

1 
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where the weight Wi is given by 



Wi= I f{dB{x,Z))hB{dB{x,Z),rB{x,Z)) ^ dx. 
J A, 

For B = B'^we have /ib2 (t, r) = 27r(t + r) (see Remark|23]l, hence if f{t) = e"^l{t < e} then 



Wi = / — ; dx. 

2vre JA,n(C,+ei?2) (x, Ci) + 

If we let e — 0, the weights converge to Wi = /j/(27rrj) if r j > and to Wi = 1 if ri = li = 0. Then the 
estimator becomes 

with li/ri interpreted as 27r if Vi = li = 0. Notice also that the outer sampling window Ty(^) then 
shrinks to W , so that in the limit only information in W is needed. The estimator Go was discussed by 
Hall |3, Chapter 5.6] (more generally, he considered estimators of EA(i2), for a given function A; Go 
corresponds to the case A = Ic). 

For B = [0, u] (with u G {±ei, ±62}), assuming (in the linear case) that G({0}) = and hence 
Ti > 0, and again choosing f{t) = £~^l{t < e}, we get hB{t, r) = 2r and 

Wi = - — / dx. 



2eri J A,r^{c,+e[Q~u]) 
This yields an estimator „ in the limit e — )• 0, which is given by 

Here, li{u) is the length of the projection of the visible part of Ci in direction u (projected onto the line 
orthogonal to u). The estimator can be improved by combining u = ei, — ei, 62, —62, 

'^l '■= ^ {^ii^ei + G/ _ei + G^^ej + Gi^^e2^ 

For applications, it would be natural to choose e = 1 which yields weights 

Wi — — , Z — i, . . . , K, 

2r,; 



and gives the estimator 



k I , I \ ~^ k 



Hence, in this case and with u = ei, information in [0, 2] x [0, 1] would be required and the estimation 
is based on the areas of the regions Ai C [0, 1]^. Of course, the estimation can be again improved by 
combining the estimators for n = ei, — ei, 62, —62 which are available if Z is observed in [—1, 2]^. 

If we do not have information from outside W, then we may use a minus sampling approach and 
replace W by the eroded window Wq^ := {x ^ W : x + eB C W}, i.e. we consider the following 
estimator 
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Another possibility would be to use the naive approach which ignores edge effects. Then we have the 
uncorrected estimator 

^ fr<\ _ W,u(C) 

where 

Jw hB[dB{x,Z r\W),rB{x,Z r\W)) 

\f B = [0, u\, then it can happen that (x, Z n W) = oo. In that case we use our convention concerning 
/ • h'j^, i.e. the points x satisfying dsix, Z n W) = oo do not contribute to rjw,uiC). Besides minus 
sampling there exist more sophisticated methods of edge correction in the statistics of spatial point pro- 
cesses. We adopt the idea of local minus sampling that was originally applied in lH to the estimation 
of the nearest neighbour distance distribution function for stationary point processes (see also [5 |). We 
use only points that are closer to Z than to the boundary of the window W. This gives the Hanisch type 
estimator 

where 

W,h(C):= / l{TB{x,Z)eC]l{dB{x,Z)<dB{x,dW)] . / ^'^''If' ^,. dx. 
Jw hB[dB{x,Z),rB{x,Z)) 

Note that for B = [0,u\ the estimators Gh and Gu coincide. 

In practical applications one has to replace in (11.51 ) the integration with respect to Lebesgue measure 
by an integration with respect to a discrete measure. This still gives a ratio-unbiased estimator of G. 

We compare the performance of the different estimators discussed above through computer simula- 
tions. We simulate a stationary planar Boolean model with spherical grains, given by (II. It . The obser- 
vation window W is the unit square [0, 1]^. The distribution G is assumed to be uniform on (0.05, 0.1). 
We approximate the integrals over W by Riemannian sums over a rectangular grid of points n W, 
where 

Lh := { {{k - l/2)h, {I - l/2)h) : A;, / G N} . 

For our purposes, we choose h = 1 /300. 

We take f{t) = e~^l{t < e} for different choices of e and compare the estimator G, given by 
(11.81) . with the estimators Go (for spherical B) and G/ (for linear B) corresponding to the limiting case 
e — )• 0. The estimators Gqe, Gu and Gh are also evaluated. For linear B = [0, u] we always combine 
the corresponding estimators for u = ei, — ei, 62, —62, this leads to a noticeable improvement. 

The radius distribution G is uniquely determined by the distribution function G{t) = G([0, t]),t > 0. 
We measure the quality of the estimators by the Kolmogorov-Smimov distance 

dKs{G,G) :=sup|G(s)-G(s)| 

s>0 



and the Cramer-von Mises distance 



/•°-^ - 9 ds 

dcvM(G,G):= / {G{s)-G{s)Y--. 



We have generated 100 independent realizations of the Boolean model Z with chosen intensity 7. 
For each realization we have determined several estimators under study. The sample means of corre- 
sponding Kolmogorov-Smimov and Cramer-von Mises distances over 100 simulations are demonstrated 
in Table [U for 7 = 25 and in Table |2] for 7 = 100. The results show that smaller values of e are more 
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^KS 




1000 • 


dcvM 


estimator 


spherical B 


linear B 


spherical B 


linear B 


G, e = 1 


0.178 


0.147 


7.921 


5.139 


G, e = 0.05 


0.172 


0.170 


7.317 


7.101 


G, e = 0.01 


0.172 


0.172 


7.295 


7.292 


Go or G; 


0.171 


0.172 


7.243 


7.257 


Gqs, £ = 0.05 


0.191 


0.177 


9.243 


7.753 


Ge„ e = 0.01 


0.176 


0.173 


7.674 


7.435 


Gu, e = 1 


0.182 


0.179 


8.389 


7.890 


Gu, e = 0.05 


0.173 


0.169 


7.480 


7.553 


Gu, £ = 0.01 


0.173 


0.168 


7.322 


7.472 


Gh, £ = 1 


0.187 


0.179 


9.003 


7.890 


Gh, £ = 0.05 


0.179 


0.169 


8.023 


7.553 


Gh, £ = 0.01 


0.174 


0.168 


7.462 


7.472 



Table 1: Sample means of distances between distribution functions computed from 100 realizations of a 
Boolean model with intensity 7 = 25 and uniform radius distribution on (0.05, 0.1). 

preferable. The limiting estimators Go and G; produced the smallest error. They are outperformed only 
in the case of smaller intensity and linear B where our estimator, given by (ll.81 l. with larger £, gives 
better results. However, this estimator uses also information from outside W. Simulation studies for 
exponentially distributed radii (not presented) show very similar results. A change of resolution h has 
only a minor influence on the quality of the estimators. For intensity 7 » 100 the deviation from the 
radius distribution increases which is intuitively clear because many balls are covered so that their radii 
are not available for the estimators. 
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